Thermodynamic phase diagram of Fe(Se .5Te .5) single crystals up to 28 Tesla 
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We report on specific heat (C p ), transport, Hall probe and penetration depth measurements 
performed on Fe(Seo.5Teo.s) single crystals (T c ~ 14 K). The thermodynamic upper critical field 
H C 2 lines has been deduced from C v measurements up to 28 T for both H\\c and H\\ab, and compared 
to the lines deduced from transport measurements (up to 55 T in pulsed magnetic fields). We show 
that this thermodynamic H C 2 line presents a very strong downward curvature for T — > T c which 
is not visible in transport measurements. This temperature dependence associated to an upward 
curvature of the field dependence of the Sommerfeld coefficient confirm that H C 2 is limited by 
paramagnetic effects. Surprisingly this paramagnetic limit is visible here up to T/T c ~ 0.99 (for 
H||a6) which is the consequence of a very small value of the coherence length £ c (0) ~ 4A (and 
£ab(0) ~ 15A), confirming the strong renormalisation of the effective mass (as compared to DMFT 
calculations) previously observed in ARPES measurements [Phys. Rev. Lett. 104, 097002 (2010)]. 
H c i measurements lead to \ a b(0) = 430 ± 50 nm and A c (0) = 1600 ± 200 nm and the corresponding 
anisotropy is approximatively temperature independent (~ 4), being close to the anisotropy of H C 2 
for T — > T c . The temperature dependence of both A (oc T 2 ) and the electronic contribution to the 
specific heat confirm the non conventional coupling mechanism in this system. 

PACS numbers: 74.60.Ec, 74.60. Ge 



I. INTRODUCTION 

The discovery of superconductivity up to 55K in iron- 
based systems [1] has generated tremendous interest. 
Among those, iron selenium (FeSei_^) [2] has been re- 
ported to be superconducting with a critical temperature 
of 8 K at ambient pressure, rising to 34-37 K under 7-15 
GPa [3] . On the other hand, the substitution of tellurium 
on the selenium site in Fei + 5(Te 3; Sei_ 2: ) increases T c to 
a maximum on the order of 14-15 K at ambient pres- 
sure [U [S] (for x ~ 0.5). This binary compound is very 
interesting as it shares the most salient characteristics 
of iron based systems (square-planar lattice of Fe with 
tetrahedral coordination) but has the simplest crystallo- 
graphic structure among Fe-based superconductors (no 
charge reservoir [6], so-called 11-structure). Moreover, 
even though the endpoint Fei+^Te [7] compound displays 
antiferromagnetic ordering, a magnetic resonance similar 
to that observed in other parent compounds (with a (1/2, 
1/2) nesting vector connecting the T and M points of the 
Fermi surface) is recovered for intermediate Te contents 
[51 |H] suggesting a common mechanism for superconduc- 
tivity in all iron based superconductors. However, in con- 
trast to iron pnictides which show only weak to moderate 
correlations, recent ARPES measurements suggested the 
existence of very large mass renormaliszation factors (up 
to ~ 20 as compared to DMFT calculations) [10] indicat- 
ing that Fe(Se,Te) is a strongly correlated metal differing 
significantly from iron pnictides. 

In order to shed light on superconductivity in these 
systems, it is of fundamental importance to obtain a pre- 



cise determination of both upper and lower critical fields 
and their anisotropy. Up to now H c2 has mainly been 
deduced from transport measurements [TTHTB] , and more 
recently by specific heat up to 14 T [J3] . As in other pnic- 
tides (see [H] and references therein), high H C 2(0) val- 
ues have been reported but, in the case of Fe(Te a; Sei_ :c ), 
strong deviations from the standard Werthamer-Helfand- 
Hohenberg model for H C 2(T) have been reported. Those 
deviations have been associated to paramagnetic limita- 
tions (so-called Pauli limit) [TIT - TIB"] . However, in presence 
of strong thermal fluctuations (see discussion below), the 
determination of H C 2 from transport measurement be- 
comes very hazardous and a thorough analysis was hence 
lacking of an unambiguous determination of H c2 from 
specific heat measurements. 

We show that the H c2 lines actually display a very 
strong downwards curvature close to T c corresponding to 
fiodH C 2/dT values rising up to ~ 12 T/K for H\\c and 
even ~ 45 T/K for H\\ab. This strong curvature, not 
visible in transport data, shows that H C 2 remains limited 
by paramagnetic effects up to temperatures very close 
to T c (up to T/T c ~ 0.99 for H\\ab). The correspond- 
ing Pauli field H p is slightly anisotropic (Hp ab /Hp C ~ 
0.8) whereas the orbital limit (H o (0)) presents a much 
stronger anisotropy 7J o (0)" ah /7J o (0)" c ^3 — 4. The huge 
H H o {<S) values (~ 130 ± 20 T for H\\c and - 400 ± 50 T 
for H\\ab) correspond to very small coherence length val- 
ues (f o &(0) ~ 15 ± I A and £ c (0) ~ 4±lA) confirming the 
large value of the effective mass previously observed by 
ARPES [10] and hence supporting the presence of strong 
electronic correlations in this system. 
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FIG. 1: AC specific heat measurements C p /T 2 as a func- 
tion of T of a Fe(Seo.5Teo.5) single crystal (sample A4) for 
IXoH = 0, 8, 12, 16, 20, 24 and 28 T (from right to left) for 
H\\ab (A) and H\\c (B). The data have been renormalized 
taking C P (T = 20K) = 3.8 J/molK. The H c2 line is deduced 
from the midpoint of the specific heat jump after subtrac- 
tion of a smooth polynomial background. (C) : Specific heat 
from relaxation data (H\\c) for the indicated magnetic fields 
(sample A5). Inset : temperature dependence of the elec- 
tronic contribution to the specific heat Ce — C p — /3T 3 — ST 5 
(solid symbols) where the phonon contribution (/3T 3 + ST 5 ) 
has been subtracted from the normal state data (see thin line 
in Fig.lC). The BCS behavior for 2S/kT c = 3.5 (solid line) 
and 2A/kT c = 5 (dotted line) are displayed for comparison. 



In addition, preliminary H c \ measurements led to con- 
tradictory results. On the one hand, Yadav et al. [TS] 
reported on rather high H c i values ~ 100 G and ~ 400 G 
for H\\c and H\\ab respectively with H c \ lines showing a 
clear upward curvature at low temperature. On the other 
hand, Bcndclc et al. |17j obtained much smaller values 
(~ 20 G and ~ 45 G for H\\c and H\\ab respectively) as- 
sociated with a clear saturation of the H C \(T) lines at low 
temperature. Finally, Kim et al. [18] reported on strong 
deviations of the temperature dependance of the super- 
fluid density (p s oc 1/A 2 oc H c \) from the standard be- 



havior, attributed to a clear signature of multigap super- 
conductivity. We present here detailed first penetration 
field measurements performed with Hall sensor arrays in 
a variety of single crystals showing very different aspect 
ratios. We hence obtained /^o-ffJi(O) = 78 ± 5 G and 
H O H l }f(0) = 23 ±3 G. The H cl lines clearly flatten off at 
low temperature but do not show the pronounced deep 
previously obtained in Tunnel Diode Oscillator (TDO) 
measurements |18j . Our TDO measurements however led 
to a similar deep which is probably due to an overestima- 
tion of the absolute AA(T) value related to spurious edge 
effects. Finally, we obtained a temperature independent 
T Hcl = H^/H^ b values ~ 3.3 ± 0.5 which corresponds 
to T\ = \ c /\ a b ~ 4.0 ± 0.8 (see below), being close to 
the 7jj o2 value obtained for T — > T c (i.e. ~ Ho ab /Ho C ). 

Finally, we confirm that A oc T 2 in both crystallo- 
graphic directions and show that the temperature de- 
pendence of C p strongly deviates from the standard BCS 
weak coupling behavior confirming the non conventional 
coupling mechanism of this system. However, the ampli- 
tude of the specific heat jump is much larger than those 
previously reported in other Fe(Se,Te) samples and hence 
does not follow the AC P vs Tj? scaling law reported in iron 
based systems [T(|l [5D] . 



II. SAMPLE PREPARATION AND 
EXPERIMENTS 

We present here specific heat, transport, Hall probe 
and Tunnel Diode Oscillator (penetration depth) mea- 
surements performed in Fei+^Seo.soTeo.s) single crystals 
grown by two different techniques. Samples A have been 
grown using the sealed quartz tube method. The samples 
were prepared from very pure iron and tellurium pieces 
and selenium shots in a 1:0.5:0.5 ratio, loaded together in 
a quartz tube which has been sealed under vacuum. The 
elements were heated slowly (100 ° C/h) at 500 ° C for 10 
h, then melted at 1000 ° C for 20h, cooled slowly down 
to 350 ° C at 5 ° C /h, and finally cooled faster by switch- 
ing off the furnace. Single crystals were extracted me- 
chanically from the resulting ball, the crystals being easy 
cleaved perpendicular to their c crystallographic axis. 
The refined lattice parameters of the Fei+^Seo.sTeo.s) 
tetragonal main phase, a = 3.7992(7) A and c = 6.033(2) 
A, are in agreement with the literature [H [17]. The real 
composition of the crystals checked by x-ray energy dis- 
persive micro- analysis using a scanning electron micro- 
scope was found to be Fei. 5(2)(Te .55(2)Se .45(2))- The 
temperature dependence of the resistivity shows a metal- 
lic behavior at low temperature as expected for this low 
level (8 = 0.05) of interstitial iron [2"T] . 

Samples of batch B were grown with the Bridgman 
technique using a double wall quartz ampoule. The in- 
side tube had a tipped bottom with a 30° angle and an 
open top. The inside wall of the outer ampoule was car- 
bon coated to achieve the lowest possible oxygen par- 



3 



TABLE I: average thickness d, width w length I or mass m of 
the samples and measurement techniques [C p = specific heat, 
R = transport, HP = Hall probe and TDO = tunnel diode 
oscillator] 



Sample 


d ((im) 


w (^m) 


I (/im) 


measured by 


Al 


50 


180 


220 


AC-C P , HP, TDO 


A2 


60 


300 


750 


HP, TDO 


A3 


65 


400 


600 


HP, TDO 


A3' 


40 


400 


300 


HP, TDO 


A3" 


40 


100 


100 


HP, TDO 


A4 




m ~ 50 /ig 




AC _ Cp 


A5 




m ~ 1.1 mg 




DC-Cp, R, 


Bl 




m ~ 0.7 mg 




DC-C P , R 



tial pressure during the growth. The Bridgman ampoule 
was inserted in a three zone gradient furnace (1000°/840° 
/700°) and lowered at a speed of 3 mm/h. At the end 
of the growth, temperature was lowered to room tem- 
perature at 50 ° C/h. Further characterizations of these 
crystals have been published elsewhere [T3]. The differ- 
ent single crystals used in this study have been listed in 
Table 1. 

The C p measurements have been performed in mag- 
netic fields up to 28 T using both an AC high sensitivity 
technique and a conventional relaxation technique. For 
AC measurements, heat was supplied to the sample by 
a light emitting diode via an optical fiber and the corre- 
sponding temperature oscillations were recorded with a 
thermocouple (sample Al and A4). In parallel, the spe- 
cific heat (sample A5 and Bl) were carried out in a minia- 
turized high-resolution micro-calorimeter using the long- 
relaxation technique. The chip resistance used as both 
thermometer and heater as well as the thermal conduc- 
tance of its leads have been carefully calibrated up to 
28T using a capacitance thermometer. Each relaxation 
provides about 1000 data points over a temperature inter- 
val of about 80% above the base temperature which has 
been varied between 1.8 and 20 K. Data can be recorded 
during heating and cooling. The merging of the upward 
and downward relaxation data provides a highly reliable 
check of the accuracy of this method. 

Electrical transport measurements have been per- 
formed on sample Bl in static magnetic fields up to 28T 
and pulsed magnetic fields up to 55T and are described 
in detail elsewhere [13] ■ We have also measured the re- 
sistivity of sample A5 with a commercial device (PPMS) 
up to 9T. 

The real part of the AC transmittivity, T' H , of sam- 
ples Al to A3" has been measured by centering these 
on a miniature GaAs-based quantum well Hall Sensor 
(of dimension 8x8 /im 2 ). The sensor is used to record 
the time-varying component B ac of the local magnetic 
induction as the sample is exposed to an ac field ~ 1 
Oe ( uj ~ 210 Hz). T' H is then defined as : T' H = 
[B ac {T) - B ac (A.2K)]/[B ac {T » T c ) - B QC (4.2K)]. The 
remanent local DC field (B rem (H a )) in the sample has 
been measured after applying a magnetic field H a and 
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FIG. 2: (color on line) (A) H — T phase diagram (batch 
A) in both H\\ab (open symbols) and H\\c (closed symbols) 
displaying the H C 2 line deduced from specific heat measure- 
ments (see Fig.l)(squares), the field corresponding to the on- 
set of the diamagnetic screening (circles) as well as the fields 
corresponding to zero resistance (downward triangles) and 
R/Rn = 0.5 (upward triangle, Rn being the normal state 
resistance). (B) same as in (A) for sample Bl. The transport 
data are taken from [13]. See Fig. 3 for a direct comparison 
between the H C 2 lines in each batch. 



sweeping the field back to zero. In the Meissner state, no 
vortices penetrate the sample and B rem remains equal 
to zero up to H a = Hf (the first penetration field). A 
finite remanent field is then obtained for field amplitudes 
larger than Hf as vortices remain pinned in the sample. 

Finally, the London magnetic penetration depth in the 
Meissner state , A, has been measured on the same sam- 
ples with a LC oscillating circuit (14MHz) driven by a 
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Tunnel Diode (TDO). The samples have been glued at 
the bottom of a sapphire rod which were introduced in 
a coil of inductance L. The variation of the penetra- 
tion depth induce a change in L and hence a shift of 
the resonant frequency Sf(T) = f(T) - f(T min ). 5f(T), 
renormalised to the frequency shift corresponding to the 
extraction of the sample from the coil A/ is then equal to 
the magnetic susceptibility. At low temperatures (typi- 
cally for T < 12K), A << d (d being the lowest dimension 

of the sample, here the thickness), and we have S &fJ = ^ 

where A is an effective penetration depth depending on 
the field orientation and R an effective dimension of the 
sample. When the magnetic field is applied along the 
c-axis, only the in-plane supercurrents are probed and 
A = \ a b, whereas A = A a b + ^A c for H//ab (w being 
the width of the sample). The effective dimension R is 
calculated following [22] . 



III. UPPER CRITICAL FIELD 

Fig.l displays typical AC measurements for both H\\c 
and H \\ab (sample A4). As shown, a well defined specific 
heat jump is obtained at T c for H = (~ 20% of the to- 
tal C p ) and this peak progressively shifts towards lower 
temperature as the magnetic field is increased (here up 
to 28 T). The H c2 line has been deduced from the mid- 
point of the Cp/T anomaly after subtraction of a smooth 
polynomial background from the raw data. As shown in 
Fig.2A, the corresponding H c2 lines present a very strong 
downward curvature for T — > T c which was not revealed 
by previous transport measurements (the same behav- 
ior is observed in all measured samples, see for instance 
Fig.2B and Fig.3A for a comparison between samples A4 
and Bl). Note that a very similar curvature has been 
reported reported very recently from C p measurements 
up to 14T [14] 

Such a curvature is a strong indication for paramag- 
netic effects and we have hence fitted the experimental 
data using a weak coupling BCS clean limit model includ- 
ing both orbital and Pauli limitations [23]. This model 
only requires two fitting parameters (plus T c ) : the ini- 
tial slope dH C 2/ 'dT\T=T a and the zero temperature Pauli 
limit H p . The results are shown in Fig.3A for sample 
A4 and Bl. As shown, very good fits can be obtained 
in both samples using very similar fitting parameters : 
HodH c2 /dT\ T=Tc - 38 ± 3 T/K and ~ 13 ± 2 T/K for 
H\\ab and H\\c respectively and n^Hp ~ 45 ± 2 T and 
~ 54 ± 4 T/K for H\\ab and H\\c, respectively. 

As previously observed in layered systems (see [24] and 
discussion in 1 
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ZfJ, 6 is actually very close to a (1 — 
T/Tc) 1 / 2 law. Strikingly, this simple behaviour is valid 
up to T/T c ~ 0.99 in our system (see Fig.3B). Such a 
dependence can be directly inferred from a Ginzburg- 
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FIG. 3: (color online) (A) comparison between the Hca(T) 
values deduced from specific heat measurements in sample 
A4 (blue squares) and sample Bl (red circles) for both H\\c 
(closed symbols) and H\\ab (open symbols) and theoretical 
values for clean weakly coupled BCS superconductors includ- 
ing both orbital and Pauli limitations (solid lines). (B) H C 2 
vs (1 — t) (sample A4) in a log-log scale (t = T/T c ) showing 
that the linear dependence of the H C 2 line rapidly crosses to 
a ~ (1 — t) ' 5 dependence. The solid lines are fits to Eq.(l) 
(see text for details). The linear slopes close to T c (~ 45 T/K 
and ~ 12 T/K) extrapolate to very high orbital limits. On 
the contrary the irreversibility line (blue circles) displays the 
(1 — t) 2 dependence characteristic of vortex melting. 



Landau (GL) expansion which leads to 
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(where H Q is the orbital field) i.e. H c2 ~ H p (^ — t) ' 5 for 
H » H 2 p /H . A fit to Eq.(l) (solid line in Fig.3B) leads 
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to fj, H p {ab - 65 T and fj, H p lc ~ 75 T, ^oH l } ab ~ 650 T 
and noH l } c ~ 170 T (sample A4) [H]. We hence have 
HoHp/H ~ 6 T for H\\ab, field which is reached for 
T/T c ~ 0.99. Fe(Se,Te) is thus a rare example of su- 
perconductor for which the upper critical field is dom- 
inated by paramagnetic effects on almost the totality 
of the phase diagram (for H\\ab). A shown in Fig.3B, 
a linear dependence is recovered very close to T c with 
[i dHt b /dT ~ 45 T/K and ^ dH o \\c/dT ~ 12 T/K, in 
good agreement with a values deduced from the BCS fit- 
ting procedure [27] . 

Those extremely high iJ Q values are related to very 
small values of the coherence lengths ^ a b — $o/27r[0.7 x 

Moffo] ~ 15 ± lA and £ c = U x {H l j c / H \\ab) ~ 4 ± 
lA which confirm the very strong renormalization of the 
Fermi velocity observed in ARPES measurements [TU] 
(see also theoretical calculations in [28 ). Indeed, one 
gets vp,ab = 7rA£ ai) //j ~ 1.4 x 10 4 m/s (A being the 
superconducting gap ~ 2 meV EO]) i-e. tivp^ab ~ 0.09 
eVA in perfect agreement with ARPES data which also 
led to hvp ~ 0.09 eVA for the a.3 hole pocket centered 
on the r point [note that the H C 2 line will be dominated 
by the band having the larger critical field i.e. the lower 
Fermi velocity] . Our measurements do hence confirm the 
strong correlation effects previously suggested by ARPES 
measurements [10] . 

An estimate of the paramagnetic field in the weak cou- 
pling limit is given by the Clogston-Chandrasekhar for- 
mula : HoHp — 2A/v / 2.g/i_B ~ 26 T in our sample (taking 
g = 2) i.e. well below the experimental suggesting that 
g ~ 1.0 — 1.2. however, it is important to note that H p 
may be increased by strong coupling effects [31] and a 
fit to the data can be obtained introducing an electron- 
phonon coupling constant A ~ 0.6 — 0.7 and g ~ 2 (still 
having an anisotropy on the order of 1.2 between the two 
main crystallographic axis) . Even if it is difficult to con- 
clude on the exact value of g, our data clearly indicate 
a small anisotropy of this coefficient (~ 1.2) supporting 
the possibility of a crossing of the H C 2 lines at low tem- 
perature. Note that this anisotropy is much lower than 
the one inferred from transport measurements (~ 4 [13] ) 
confirming that the large apparent anisotropy of g de- 
duced from those measurements is an artifact, probably 
related to the anisotropy of flux dynamics (see discus- 
sion on the irreversibiliy line below). The anisotropy 
of the upper critical field is then strongly temperature 
dependent rising from H^/H% 2 ~ Hf/H° ~ 0.8 for 
T 0, reflecting the small anisotropy of the g factor, to 
H ci/ H c2 ~ Hf/H C Q - 3.5 - 4 close to T c , reflecting the 
anisotropy of the coherence lengths (see Fig. 7). 



IV. IRREVERSIBILITY LINE 

The small £ values associated to large A values 
(A a &(0) ~ 430 nm (see below and [T7J) lead to strong 
fluctuation effects hindering any direct determination of 
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FIG. 4: Transport and AC transmitivity measurements as 
a function of T for the indicated magnetic fields (H||c) in 
Fe(Seo.5,Teo.5) single crystals. In the inset : comparison be- 
tween transport and specific heat data for fioH — and 6 
T (\\c) emphasizing that the midpoint of the specific heat 
anomaly does not correspond to any characteristic tempera- 
ture in R(T) for H £ 0. 



H c2 from either transport of susceptibility measurements. 
These fluctuations can be quantified by the Ginzburg 
number d = (k B T c /e £c) 2 /$> where e (= ($o/47rA a6 ) 2 ) 
is the line tension of the vortex matter. One hence ob- 
tains eo£ c ~ 40 K (as a comparison e £c ~ 200K in 
cuprates) and Gi ~ 10 -2 which is very similar to the 
value obtained in YBa2Cu30y_5 or NdAsFe(Oi_ x F I ) (so 
called 1 Ill-phase, see [32] and references therein) clearly 
showing that thermal fluctuations are very strong in this 
system. 

To emphasize this point, we have reported in Fig. 2, 
the temperatures corresponding to both B — > and 
R/Rn — 0.5 deduced from transport measurements up 
to 9T for sample A4 (see also Fig. 4) and even up to 50T 
for sample Bl (see [T3]) (Rn being the normal state resis- 
tance). As shown, none of those lines present the strong 
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downward curvature obtained in C p measurements. On 
the contrary, the R/Rn — lines vary almost linearly 
with T with d^HjdT ~ 11 T/K and ~ 5 T/K for H\\ab 
and H\\c, respectively in agreement with previous mea- 
surements [TTJ [T2]. However, as pointed above, these 
lines do not correspond to any thermodynamic criterion 
and discussions of the corresponding lines should hence 
be taken with great caution. Moreover whereas the mid- 
point of the specific heat coincides with the R = tem- 
perature for H = in sample A4, this midpoint rather 
lies close to the R/Rn = 0.5 point in sample Bl clearly 
showing that neither of those two transport criteria can 
be associated with the H c2 line. 

Similarly, as previously observed in high temperature 
cuprates and 1111-pnictides [52] , the onset of the dia- 
magnetic response (T' H — > 0) also lies well below the the 
H C 2 line, (see Fig.2A and Fig. 4). Indeed, this onset is 
related to the irreversibility line above which the system 
is unable to screen the applied AC field due to the free 
motion of vortices. This irreversibility line is then ex- 
pected to lie close to the R = line. As shown in Fig. 4 
the onset of diamagnetism actually differs slightly from 
the onset of resistivity. This difference is much probably 
related to different voltage-current criteria (the magnetic 
screening corresponds to much smaller electric fields but 
requires higher currents) but both lines present the pos- 
itive curvature characteristic of the onset of irreversible 
processes. Note that, as expected for vortex melting (for 
a review see |33j). the irreversibility line (here defined as 
the onset of T' H ) varies as : H irr oc (1 — T/T c ) a with 
a ~ 2 (see Fig.3B). A similar curvature has also been 
reported by Bendele al. [T7] for the irreversibility field 
deduced from magnetization measurements. 

V. LOWER CRITICAL FIELD 

The first penetration field has been measured on a se- 
ries of Fe(Seo.5Teo.5) samples with very different aspect 
ratios (see Table 1). To avoid spurious effects associ- 
ated to strong pinning preventing the vortex diffusion 
to the center of the sample [34 Hf has also been mea- 
sured on several locations of the same sample. The in- 
set of Fig. 5 displays typical examples on sample A3' (2 
positions) and A3". In samples with rectangular cross 
sections, flux lines partially penetrate into the sample 
through the sharp corners even for H a < H t but remain 
"pinned" at the sample equator. The magnetization at 
H a = Hf is then larger than H c \ and the standard "el- 
liptical" correction for H c \ (= — N) where N is 
the demagnetization factor) can not be used anymore. 
Following |35j . in presence of geometrical barriers, Hf is 
related to H c \ through : 



tanh(Wad/w) 

where a varies from 0.36 in strips to 0.67 in disks (d and 
w being the thickness and width of the sample, respec- 
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FIG. 5: Temperature dependence of the lower critical field 
(H c i) deduced from Hall probe measurements for £f||a& (open 
symbols) and H\\c (closed symbols) in sample Al (squares), 
A2 (circles), A3 (lozenges) and A3" (triangles), see Table 1 for 
sample dimensions. Solid lines are guides to the eyes, crossed 
squares are muon relaxation data from [23] and open crosses 
TDO data. In the inset : remanent field B rem as a function 
of the applied field in sample A3' (center and edge) and A3" . 



tively. To reduce the uncertainty associated with the a 
value as well as the d/w ratio in real samples of irregular 
shape, five different samples with different aspect ratios 
have been measured (see Table 1). Sample A3' has been 
cut out of sample A3 and finally A3" out of A3' in or- 
der to directly check the influence of the aspect ratio on 
Hf. The corresponding Hf values are reported in the 
inset of Fig. 7 together with the theoretical predictions 
from Eq.(2) taking /io-^ci = 78 G (the predictions for a 
an standard "elliptical" correction are also displayed for 
comparison) . 

The lower critical fields (fioH^, fXoH^) are then re- 
lated to the penetration depth (A c , X a b) through : 

HoH c cl = -^(Ln( K ) + c( K )) (3) 

47TA a bA c 

where k — \ a b/£ab, K * = A c /£ a b and c(k) is a k dependent 
function tending towards ~ 0.5 for large k values. Taking 
fi H^ 2 (0) = 0.7 x n H o - 130T, and H c cl = 78 ± 5 G one 
gets A a b(0) ~ 430±50 nm, which is in fair agreement with 
muons relaxation data [17l [36] . This very large A value 
confirms the general trend previously inferred in iron 
pnictides (see for instance |17j and references therein) 
pointing towards a linear increase of T c vs 1/ \ 2 ab as ini- 
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FIG. 6: Temperature dependence of X a b (solid symbols) and 
A c (open symbols) deduced from the frequency shift in TDO 
measurements (sample A3", Table 1). Inset : temperature 
dependence of the superfluid density p™° (T) / Ps DO (0) = 
1/(1 + AA a6 (T)/A ab (0)) 2 taking A a6 (0) = 430 nm and R = 
14/xm (solid symbols) (i.e. following [22] , see correspond- 
ing AA a (,(T) values on the main panel) or R = 70/im (open 
symbols). The average H c i (T)/H c \ (0) curve (see Fig. 5) is 
displayed as the thick solid line. 



tially proposed in cuprates by Uemura et al. [37]. For 
H\\ab, no correction was introduced (flat samples) and 
one hence obtains poH^i = poHj b = 23 ± 3 G leading to 
A c - 1600 ± 200 nm (taking p, H$ = 0.7 x p Hf ~ 460 
T) 

As shown in Fig. 5, H C \(T) clearly shows a saturation 
at low temperature. As a comparison we have reported 
on Fig. 5 the temperature dependence of the superfluid 
density deduced from muons relaxation data |36j and 
Ps DO ( T ) measurements [TS]. Both H cl and pf R (T) 
curves are similar but do not reproduce the important 
shoulder at 5 K of the superfluid density. This shoulder 
has been interpreted as a clear signature of multi-gap su- 
perconductivity and as a failure of the clean limit s-wave 
(including s±) pairing |18j . Our measurements do not 
support this interpretation. 

In order to shed light on this discrepancy, we have 
performed TDO measurements on each of the samples of 
Table 1. As described in sec. II, A c and A ab were deduced 
from the effective penetration depth A measured for both 
H\\ab and H\\c. As shown in Fig. 6 (sample A3"), both 
AA a b and AA C are proportional to T" with n close to 
2, in good agreement with previous measurements for 
H\\c [18] (the same temperature dependence has been 



obtained for all samples). The TDO data then require the 
introduction of the value of A ab (0) to convert the AA(T) 
data into p T s DO (T) / p T s DO {$) = 1/(1 + AA Qb (T)/A ab (0)) 2 . 
Introducing A ab (0) ~ 430nm and taking R ~ 14/xm 
(from [H]), Ps DO {T) shows a change of curvature around 
5K, very similar to the one previously reported in [18] (see 
inset of Fig.6). A similar discrepancy has already been 
observed in MgCNi3 and interpreted as a reduction of the 
critical temperature at the surface of the sample due to a 
modification of the carbon stoechiometry |38j . However, 
such an explanation is not expected to hold here as single 
crystals were extracted mechanically from the bulk. 

It is important to note that the temperature depen- 
dence of the superfluid density is very sensitive to the 
absolute value of AA and, although very similar to the 
one reported by Kim et al. |18j . the amplitude of 
AA Q b/T 2 ~ 40A/K 2 observed in our samples is much 
larger than the one reported recently by Sarafin et al. 
[T3] (~ lOA/K 2 ). Similar discrepancies in the absolute 
amplitude of AA have also been reported in other pnic- 
tides j3H] and have been attributed to complications from 
rough edges which may lead to an overestimation of AA. 
Dividing the absolute AA ab by a factor ~ 5 (i.e. taking 
R = 70 /jm for H\\c instead of 14 pm) actually leads to 
a very good agreement between TDO and H c \ data (see 
Fig.6) hence indicating that this value has probably been 
overestimated due to an underestimation of the effective 
dimension R in presence of rough edges. 

Very similar temperature dependences of H c \ were ob- 
tained in both directions (see Fig. 5) leading to a (almost) 
temperature independent anisotropy of H c i : Th c1 ~ 
3.4±0.5 and hence T A = A c /A ab = [H^/H^]x [{Ln( K *) + 
c(k))/(Im(k*) + c(«))] ~ T Hcl x 1.2 - 4.1 ± 0.8 (see Fig. 
7). This value is hence very close to the one obtained 
for H c2 close to T c as T Hc2 {T T c ) ~ T Ha = U/L 
(see Fig. 7). Similarly, very similar temperature depen- 
dences have been obtained for AA C and AA ab (with 
AA C ~5x AAab up to T ~ T c ) again suggesting a weak 
temperature dependence of this anisotropy. Finally, this 
value is also close to the one obtained for the irreversibil- 
ity field deduced from the onset of diamagnetic screening. 



VI. FINAL DISCUSSION 

The value of the normal state Sommerfeld coefficient 
(7jv) m Fe(Se,Te) compounds remains debated as values 
ranging from - 23 mJ/molK 2 (40] to - 39 mJ/molK 2 g] 
have been obtained. For non superconducting samples, it 
has even been shown recently [41] that rises rapidly 
for x < 0.1 reaching ~ 55 mJ/molK 2 for 0.1 < x < 
0.3. Even though our maximum field (28 T) is too low 
to fully destroy superconductivity down to K hence 
hindering any precise determination of 7^, it is worth 
noting that a 7^ value on the order of ~ 39 mJ/molK 2 
is incompatible with the entropy conservation rule in our 
sample. A reasonable fit to the data (solid line in Fig.lC) 
assuming that C p /T = j N +l3T 2 +5T A for 20 > T > 12 K 
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FIG. 7: Temperature dependence of the anisotropy of the 
upper (open circles), and lower (solid squares) critical fields 
and irreversibility line (onset of diamagnetic response, open 
crosses). The Th c1 data have been multiplied by a factor 1.2 
to display Y\ = 1.2 x Th b i- In the inset : first penetration 
field Hf as a function of the aspect ratio : 2d/ (w + I). Thick 
solid lines are theoretical predictions in presence of geometri- 
cal barriers using see Eq.(2) for disks (upper line) and strips 
(lower line); the thin line correspond to a standard "elliptical" 
correction (no surface barrier). 



and noH = 28 T leads to — 23 ±3 mJ/molK 2 in good 
agreement with the value obtained by by Tsurkan et al. 
[40] , This 7tv value is also in fair agreement with the one 
deduced from ARPES measurements (~ 30 mJ/molK 2 
|10j). Similarly, the Debye temperature (873 ~ 143 K) is 
in reasonable agreement with the one previously reported 
in both Fe(Se .67Te .23) (6 D ~ 174 K [4 ) and Fei. 05 Te 
{Q D ~ 141 K [32]). 

The electronic contribution to the specific heat 
{CJT = C p /T - j3T 2 - ST 4 ) is then displayed in the 
inset of Fig.l together with the theoretical prediction 
for a single gap BCS superconductor in the weak cou- 
pling limit (i.e. taking 2A/kT c ~ 3.5, thin solid line). 
As shown, this standard behavior largely overestimates 
the experimental data at low temperature suggesting the 
presence of a much larger gap. A reasonable agreement to 
the data is obtained assuming that 2A/kT c ~ 5 (dotted 
line). However, even though some indication for the pres- 
ence of a large gap were obtained by fitting either fiSR 
[36] or optical conductivity [43] data, the corresponding 
gap value (~ 3 meV) is much larger than the value ob- 
tained by spectroscopy (~ 1.8 — 2 meV [221 [SO])- More- 
over those former measurements also suggest the presence 



of a much smaller gap which is not present in our specific 
heat meaurements. 

Some evidence for nodes (or for deep gap minima) in 
Fe(Seo.5Teo.5) has been suggested by four fold oscilla- 
tions in the low temperature specific heat for H\\c [45] . 
However, despite the high resolution of our AC tech- 
nique and the very good quality of our samples (the 
specific heat jump at T c is slightly larger than in [4"5] ) 
we did not observe these oscillations in our samples (i.e 
AC P {9) I C p < 10~ 3 ). Nodes are also expected to show 
up in the field dependence of the Sommerfeld coefficient 
(7 (if)) which is then expected to vary as H a with a < 1 
(a = 0.5 for the so-called Volovik effect for d-wave pair- 
ing with line nodes whereas a ~ 1 for classical single gap 
BCS systems). We have hence extrapolated the C e (H) /T 
data to zero using either a BCS formula (see discussion 
above, C e /T — j(H) cx exp(— A(H)/kT) in our tempera- 
ture range) or a phenomenological second order polyno- 
mial fit. Both procedure led to a concave curvature for 
j(H) with a ~ 1.5 ± 0.3 for H\\c and a ~ 2.2 ± 0.6 for 
H\\ab. This concave behavior can be attributed to the 
effect of Pauli paramagnetism on the vortex cores [46] 
(see [171 [35] for experimental data in heavy fermions) 
hence clearly supporting the importance of these effects 
in Fe(Se .5Te .5)- 

Finally note that it has been suggested that AC P /T C 
could be proportional to T 2 in iron pnictides [HI I2U] 
due to strong pair breaking effects [49] with AC p /Tj? ~ 
0.06 mJ/molK 4 . One hence would expect an anomaly 
ACp/T c ~ 12 mJ/molK 2 at T c in our system which 
is clearly lower than the experimental value ~ 40 ± 5 
mJ/molK 2 . Similarly, it has been suggested that the ini- 
tial slope of the H C 2 line could scale as nodH^/dT ~ 
0.2 x T c (T/K) but, again, this scaling docs not hold 
in our sample for which nodH^/dT ~ 12 T/K. Finally 
note that the temperature dependence of the superfluid 
density (see discussion above) supports the AX a b/T 2 ~ 
lOA/K 2 value obtained by Serafin et al. [H] which is also 
much smaller than the one suggested from the scaling of 
[5U] : A\ ab /T 2 - 8.8 x 10 4 /T C 3 ~ 32A/K 2 . 



VII. CONCLUSION 

In summary, 

(i) Precise determinations of the H c2 lines from C p 
measurements led to a very strong downward curvature, 
similar to that observed in layered systems. 

(ii) The temperature dependence of the upper critical 
field and the field dependence of the Sommerfeld coeffi- 
cient both indicate that H C 2 is limited by strong param- 
agnetic effects with fj, H p — 45 ± 2 T and — 54 ± 4 T for 
H\\ab and H\\c, respectively. 

(iii) The very small value of the coherence length 
£af>(0) ~ 15 A confirms the strong renormalisation of the 
effective mass (compared to DMFT calculations) previ- 
ously observed in ARPES measurements [TU] and associ- 
ated strong electron correlation effects. 7/v is estimated 
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to ~ 23 ±3 mJ/molK 2 in fair agreement with the ARPES 
value. 

(iv) The anisotropy of the orbital critical field is esti- 
mated to be on the order of 4 hence leading to a £ c (0) 
value smaller than the c lattice parameter. 

(v) Neither the temperature dependence of A nor that 
of the electronic contribution to the specific heat follow 
the weak coupling BCS model (an BCS dependence with 
A/kT c ~ 5 remains possible) but no evidence for nodes 
in the gap is obtained from the field dependence of the 
Sommerfeld coefficient. We did not observe the fourfold 
oscillations of the low temperature specific heat previ- 
ously obtained by Zeng et al. |45j . 

(vi) The amplitude of the specific heat jump 
ACp/T c ~ 40 ± 5 mJ/molK 2 is much larger than that 
previously observed in Fe(Se,Te) and does not follow the 
AC P /T C 3 inferred in iron pnictides. Similarly neither the 
slope of the H C 2 line nor the absolute value of AA(T) 
obey the scaling laws previously proposed for iron pnic- 
tides ga hjj . 

(vii) A ab (0) = 430±50 nm and A c (0) = 1600±200 nm, 
confirming the very small superfluid density previously 
observed in iron pnictides. The corresponding anisotropy 



is almost temperature independent with T\ ~ T h 2 (T — > 

(viii) These large A values associated to small £ val- 
ues lead to a very small condensation energy eo£ c ~ 40K 
and hence to large fluctuation effects hindering any de- 
termination of H C 2 from either transport or susceptibil- 
ity measurements. A detailed analysis of the influence of 
these fluctutations on the specific heat anomaly will be 
presented elsewhere. 

(ix) The strong upward curvature of the irreversibil- 
ity line (defined as the onset of diamagnetic screening) : 
Hi rr oc (1 — T/T c ) 2 strongly suggests the existence of a 
vortex liquid in this system. 
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M.Konczykowski from the Laboratoire des Solides Ir- 
radies, Palaiseau for the development of the Hall sensors 
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